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In those charge-coupled devices (CCD's) which have regions under 
each electrode which are substantially free of externally applied 
tangential electric fields, charge motion takes place by space charge 
assisted diffusion, and this relatively slow process represents a limita- 
tion to the operating frequency of CCD's with large plates. In this 
paper, subject to certain approximations, the equations of motion for 
CCD charge transfer have been derived, yielding a nonlinear diffusion 
equation. The solution of this equation by a stable finite difference 
scheme is described, and the solutions are applied to predicting the 
operating characteristics of CCD's. The results in synopsis are: 

(i) The n-channel devices will have lower losses at a given frequency 

than the p-channel devices, and a higher upper frequency limit. 

(it) Higher amplitude signals (more charge) yield lower losses. 

(Hi) Losses can be reduced an order of magnitude by using zero's 

which carry a substantial amount of charge. For example, with 

2-MHz clocking of a 25-nm plate (pad), p-channel 2-phase 

device, a 2-volt one is diminished by 4 percent per transfer with 

empty zero's, but with 1-volt zero's, the diminution is 0.26 

percent per transfer. 

(iv) Reasonably efficient CCD operation should be possible up to the 

50-MHz range using contemporary design tolerances. 

(v) Diffusion is important for reaching high transfer efficiencies. 

The frequency limitations described in this paper can be overcome 

by using a structure in which the distance between the electrodes 

and transferring charge is comparable to the electrode width and 

spacing. 

i. introduction 

The charge-coupled device, as conceived by Boyle and Smith 1 and 
realized by Araelio, Tompsett, and Smith, 2 - 3 consists of a series of 
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metal oxide semiconductor (MOS) capacitors that are driven by clock 
pulses into deep depletion. In this condition, the MOS structures are 
capable of holding minority carriers in the potential wells beneath 
the plates (pads) . In order to establish a transfer of charge from one 
plate to another, it is necessary to make the potential well beneath 
the second plate deeper than that beneath the first. This is illustrated 
in Fig. 1. In the absence of charge, the potential configuration beneath 
any one plate would be essentially flat, and finite electric fields would 
exist only in that span between two CCD electrodes. When charge is 
present, its transfer is driven predominately by the electrostatic forces 
associated with the presence of the charge in the CCD and by the 
thermal forces responsible for diffusion. It is the purpose of this paper 
to describe the transport of charge under the influence of these forces 
and, using this description, to make predictions concerning the opera- 
tion of CCD's. In order to do this, it will be necessary first to derive 
the equation governing the transport of charge and solve it. Then the 
solutions will be applied to some particular CCD situations. 

II. DERIVATION OF TRANSPORT EQUATION 

The derivation of the transport equation will utilize a number of 
approximations. The first of these is a linear approximation of 
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Fig. 1— A diagram of a CCD cross section and surface potential diagrams 
illustrating charge storage and charge transfer in CCD's. 
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the surface potential versus applied voltage relationship. The one- 
dimensional relationship between surface potential f and applied 
voltage V A is described by this equation: 



= V + V - VV 2 + 2V V' 



<7 

c 

(1) 



V = 8 2 eN D -^ 



^ ox - 



5 

In this equation Vfb is the MOS flatband voltage, and q is the surface 
density of mobile charge stored at the interface between the silicon 
and silicon dioxide. The constant V n depends upon the oxide thickness 
8, the doping density N D , the dielectric constant e 0J of the oxide, and 
the dielectric constant e„ of the silicon. In the subseqeunt calculations, 
however, an approximation to this relationship will be used which is 
linear in V A and q. This approximation is 

<P = <Po + i(v A - 7™ + ^-)- (2) 

In this case y is a number, typically between 0.7 and 1, which matches 
equation (2) to equation (1) over the range of anticipated operation. 
Such a match is shown in the example of Fig. 2, which shows the sur- 
face potential versus applied voltage reduced by the flatband voltage 
and q/C or , both in the precise form and in the linear approximation. 
It may be seen in this case that the linear approximation is rather 
good over the range from approximately 4 to 15 volts. If the electric 
field parallel to the Si-SiOa interface is calculated from the approxi- 
mate equation, it is found to depend upon the derivative of q. There 
is, however, another term to the tangential electric field; this term 
arises from the electrostatic repulsion of the carriers. In a two- 
dimensional approximation, the tangential electric field at x will have a 
component due to the summation of the fields of line charges located 
at other points .r / along the surface. Because there is a metallic elec- 
trode on the surface of the Si0 2 , each elemental line charge has an 
image, and the result is a shielding of the tangential field. Taking 
account of the two different dielectrics and the image charge, the 
repulsion field E K can be written as this integral: 
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Fig. 2 — The exact variation of deep depletion surface potential in an MOS with 
1000A of Si0 2 over a substrate with 5 X 10 14 ionized impurities per cm 3 , and 
the linear approximation * cc yV'. 



EM = *<€., + e.) / g^')^ " {x JV}* + 4^} 



dx' 



(3) 



Since the principal contribution to this integral comes from the small 
region within a few 8 of x, it is appropriate to expand q(x') as a 
Taylor series about x. If this is done, and the limits of integration 
are extended over all x', all even terms in the expansion vanish by- 
symmetry, and to the extent that d*q/dx* and higher derivatives are 
negligible, the repulsion field becomes a local function of dq/dx. 



E R (x) = - 



25 68 



e ox + e, dx 



(4) 



If one takes the gradient of equation (2) and adds the repulsion field 
in equation (4) , the total electric field acting on the charge under the 
CCD plate is given by 



E T = -S 



dq 
dx 



(5) 
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Typical values of both y and the net elastance S are presented in Figs. 
3 and 4 as functions of the oxide thickness 8. Ordinarily 8 is approxi- 
mately equivalent to the reciprocal of the oxide capacitance; conse- 
quently, Fig. 4 has been plotted in such a way that the value of S is 
compared with the oxide capacitance. 

The transport equation for CCD's will be based on the divergence 
equation for current 



using the relationship 



^= -V.J 

at v J > 



dx 



(6) 



(7) 



for the value of the current density J. In this equation D is the 
diffusion coefficient appropriate to the surface, and n is the surface 
mobility. Combining equations (5), (6), and (7) gives the basic equa- 
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Fig- 3 — The proportionality constant 7 as a function of oxide thickness, using 
doping density as a parameter. 
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Fig. 4— The effective surface capacitance, 1/8 (the reciprocal of elastance S) 
normalized by C oz and presented as a function of oxide thickness. 



tion for charge transport in a CCD : 

at-^toJ + a* a? + * a* 8 

This may be simplified somewhat by noting the Einstein relationship 
between diffusion and mobility. 



(8) 



n kT 



(9) 



Here kT 7 and e all assume their traditional values as Boltzmann's 
constant, absolute temperature, and electronic charge. Using this sub- 
stitution, one reaches this equation: 



^-nfli?;) - l /<•-/■;.:'-''' a- 



(10) 



dt r ""\da;/ ' r " 2 to 2 e M ax 2 " 
The last term on the right of equation (10) is simply the diffusion 
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term. The resemblance between this diffusion term and the other term 
involving the second partial of charge is quite apparent, because the 
product Sq, like nT/e, is a voltage. 

Without boundary conditions, the problem is only half specified. 
The physical situation on which the boundary conditions are based is 
shown in Fig. 5. In Fig. 5a, a CCD plate is shown storing charge. The 
potential beneath the plate and the charge density are both uniform. 
In Fig. 5b, another potential well has been impressed immediately 
adjacent to the first potential well. The second well is considerably 
deeper than the first; as a consequence, charge will start to flow from 
the first potential well into the second under the control of equation 
(8). No charge, however, flows in the opposite direction. Between the 
two potential wells there is an abrupt step in potential. This means 
that there is an extremely high electric field, and charge in that 
vicinity will move with a very high velocity. Using the preceding facts, 
it is possible to approximate the physical boundary conditions in the 
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Fig. 5 — An illustration of the length definition and the surface potential con- 
figuration in equilibrium, and at I — 0, when the transfer is initiated. 



1728 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUGUST 1971 

following way. 

te-* " x ^° (id 

q = at x = L. 
Both of these boundary conditions apply for all time. The first ex- 
presses the fact that current is unable to flow to the left in the device, 
and the second presumes either an ideal sink for charge at x = L, or 
alternatively that charge moves away from that point at an infinite 
velocity. The initial condition is a value of q at time t = for all 
points except x = L; it will be assumed that the charge is evenly 
distributed with the magnitude q . 

Before proceeding with the solution to this equation, it will be con- 
venient to scale the variables. The scaling can be effected by applying 
the following definitions: 

r = 4" = 1 (12a) 

L Tq 

MTo 
W = 1 VOlt 



y = 



I (12b) 



u = & (12c) 

W 

Time is now represented by the variable t, which scales time against 
t , and y is the dimensionless variable representing length. In the 
definitions of t and t , the voltage unit Uq is introduced. If this equa- 
tion had a natural voltage unit, Wo would be defined as that unit. 
However, there is no natural voltage scale to equation (10), even 
though the thermal voltage K T/e appears in one term ; consequently, 
it has proven convenient to define Wo to be unity. The quantity u 
representing charge comes about because the product of the elastances 
and charge q is equivalent to a voltage. In the solutions which follow 
t will typically range from lO" 2 to 10 2 , y will range from to 1, and u 
will range from 10 -3 to 10. 

III. SOLUTION OF THE EQUATION 

The scaled equation has the form 

u r = ul + (u + a)u vu (13) 
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where a = xT/euo is a constant. The problem may be simplified by 
letting w = u + a to obtain 

w r = w 2 u + WyyW (14) 

^ y ^ 1 
g r ^ oo 



subject to 



W,(0, r) = 
w(l, t) = a 

»fe,0)-l U < + ' 0S " <1 (15) 

k y = l 

where w { is determined by the initial uniform charge density q . 

The nonlinear parabolic initial boundary value problem given by 
equation (14) and equation (15) may be solved using a finite differ- 
ence scheme. The scheme comes from a more general technique for 
solving systems of coupled nonlinear parabolic equations which results 
in linear difference equations. This more general scheme will be pub- 
lished at a later date. 

Let h be the space mesh size and At be the time step. Let w) = 
w(jh, kAr). The indices j and k describe position and time respectively. 
The difference equation is then: 

«?" - w) 

At" 

" 2 L 2h 2h + Wi h 2 

<, - wU will ~ <-\ | » w)l\ + tgKj - 2<' 1 (m 
+ 2h 2h + *' h 2 J ' (16) 

for j = 1, • • • , J — 1, where h = l/J . The boundary conditions reduce 
to w k + l = w>* +1 and w k f x = a. 

This is a direct generalization of the Crank-Nicholson scheme for 
solving the heat equation, w t = w xx . For a description of that scheme 
and its properties, see Ref. 4, pp. 185-193. 

The most important properties of (16) come from the time and 
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space averaging which is the hallmark of the Crank-Nicholson scheme. 
The left-hand side of (16) is close to the value of w t (jh, {k + 1/2) At) 
and the right-hand side of (16) has been averaged to be close to 

(w u (jh, (k + 1/2)M)) 2 + w(jh, (k + l/2)M)w M (jh, (k + 1/2) A*). 

The averaging has been done in such a way as to make w)\\ , w)t\ 
and w) +1 appear linearly. This averaging gives (16) two nice properties. 
First, if a solution of (14) is inserted in (16), and everything is expanded 
in a Taylor series about (jh, (k + 1/2) M), then equality of the right- 
and left-hand sides follows up to terms of order 0(h 2 ) and 0(At 2 ). That 
is, the scheme is accurate to second-order in both time and space. 
In practice, this means that the difference between the computed 
w) and w(jh, kAt), the true solution at the mesh points, will be 0(h 2 ) 
+ 0(At 2 ). Second, (16) is a linear system of equations for the to} +1 . 
In fact, equation (16) is equivalent to 

+ <'(i - f£ («;♦. + wf., - m> + ^ wk ) 

+ w)v{^ (wU - wU) -§?«})-«}• < 17 > 

This has the form 

A,w t+1 + Btw, + O,-, = Di j = 1, 2, • • • , J - 1, (18) 

with wj = a and w<> = Wi , reflecting the boundary conditions. This is 
a tridiagonal system of equations for the Wj and may be solved quite 
efficiently using Gaussian elimination (see Ref. 4, pp. 198-201). The 
method is easily described. Let 

Wi = E iU > i+1 + F t (19) 

for j = 0, 1, • • • , J - 1. Replacing Wj-i by ff,_it* + F,_i in (18) 
gives w>j + i in terms of Wj . Comparing this relation with (19) gives 

'* B, + 0,0,-, 

for j = 0, • • • , J - 1. Since w»i = w = #0^1 + F we see that E = 1 
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and F = 0. This completely determines w , • • • , Wj using (20) for 
; = 1, ••-,</- 1, and then (19) for j = J - 1, • • • , 0. 

IV. RESULTS AND INTERPRETATION 

In applying the results of this analysis to CCD operation, the most 
important operating characteristic is the amount of charge remaining 
behind after transfer from one plate to the next in a limited time. 
The remaining charge is plotted as a function of time in Fig. 6, using 
two different initial conditions, Sq = 10 volts and Sq = 2 volts. The 
nonlinear equation which has been solved reduces essentially to the 
diffusion equation after most of the initial q has been dissipated, and 
the solutions in this domain approach straight lines of log (charge) 
versus time when the average value of Sq is much less than K T/e. 
However, at short time r <, 1 the solution is effectively driven by the 
charge so the time rate of change for Sq = 10 volts is roughly five 
times the time rate of change associated with Sq = 2 volts. Between 
these two domains is a transition region where the average value of Sq 
lies between 0.2 «T '/e and 5 kT/c Here the curves are very nearly 
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Fig. 6 — The fraction of the charge remaining under the first plate as a function 
of the normalized time r = t/r . The abscissa is plotted as t 1/2 , which spreads 
the region where t < 1, and leads to straight lines over the transition portion of 
the characteristic. 
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straight lines of log (charge) versus A The different regimes reflect 
a feature mentioned when the equation was scaled; i.e., there is no 
single natural potential to use in forming the time scale. For very 
short times, the natural potential is the initial value of Sq. For times 
much greater than one, the natural potential is kT/b. Figure 7 shows 
the transfer of charge with and without the aid of diffusion. It is clear 
that the last 0.1 volt of scaled charge receives a considerable boost 

from nT/e. 

In Fig. 6 it can be seen that in the vicinity of t = 10 the charge 
remaining diminishes to approximately 1 percent of its initial value, 
and in the vicinity of t = 30 the charge diminishes to 0.1 percent. 
These results are meaningless without some concept of the size of t. 
Figure 8 helps to alleviate this shortcoming by showing t as a func- 
tion of the plate width L for n-channel and p-channel devices. The 
n-channel mobility has been taken to be 750 and the p-channel mo- 
bility 150 cm 2 /volt-second. B From Fig. 8 it is possible to see that the 
values of t„ range from approximately 1 ns to 1 /xs for plates ranging 
from 10 to 100 jtun in width. To take a specific example, a bti-pm plate 
in a p-channel device would have a t„ of one-sixth of a microsecond. 
If operation were desired at a value of t = 10, aiming for a transfer 
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Fig. 7 — The fraction of the charge remaining for initial values of Sq n = 1 volt 
with room temperature diffusion and no diffusion (xT/e = 0). 
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Fig. 8 — A plot of the time scale factor t p as a function of the CCD plate 
length L. 



efficiency near 99 percent, then the operating frequency of a 2-phase 
CCD would be roughly 300 kHz, and the same loss per transfer 
would be achieved in a 3-phase CCD at 200 kHz. In an n-channel 
device with the same loss, those frequencies would increase to 1.5 and 
1 MHz respectively. 

In order to illustrate the effect of frequency, another example is 
shown in Fig. 9, where the loss per transfer is plotted as a function 
of clock frequency for a p-channel 2-phasc CCD operating with square 
waves. In this type of operation, higher signal levels result in higher 
transfer efficiencies. The plate width in this case was chosen to be 
25 ju,m, and at 1 MHz the loss per transfer ranges from about 0.3 
percent for Sq„ = 10 to 1.5 percent for Sq„ = 2 volts. 

The figures quoted so far must be regarded as worst-case loss 
figures. They apply to an isolated cluster of charge being propagated 
through a CCD. It will be possible to reduce this loss by approxi- 
mately an order of magnitude by arranging the information propa- 
gated through the device in such a way that no information ever 
leads to a totally empty CCD cell. This will be effective, because for 
long times, i.e., t > 10, the charge remaining behind is substantially 
independent of the amount of charge initially put under a pad. This 
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Fig. 9— The fractional loss per transfer of a P-channel 2-phase CCD, with a 
field-free length L of 25 ^m, plotted against clock frequency with the initial one 
signal as a parameter. The zero's have no charge. 

is illustrated by Fig. 10, which shows the absolute magnitude of 
remaining charge as a function of the charge which was initially 
placed in the device. In this figure, the amount of charge is described 
as its voltage equivalent Sq. The impact on device loss is illustrated 
in Fig. 11, where the loss per transfer for the same device shown in 
Fig. 8 is plotted as a function of the charge existing in a zero when a 
one contains 2 volts. It is presumed that the information carried is an 
alternating series of one's and zero's. This represents a worst-case 
condition when nonempty or "fat" zero's are used. From this figure, it 
is readily seen that a 2 to 1 one to zero ratio leads to practically an 
order of magnitude reduction in loss at a given frequency. Operation 
in this way requires the use of a threshold detector to ascertain 
whether a one or a zero was present; the success of such operation 
depends in large part upon the realization of a practical threshold 
detector. 

The solution of the nonlinear diffusion equation also gives charge 
distributions in the CCD as a function of time. The charge distribu- 
tion at selected times is shown in Fig. 12. Initially the left-hand end 
is almost totally undisturbed by the presence of an accepting well at 
the right-hand end of the CCD, but then as charge transfers out of 
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Fig. 10 — The remaining charge at several times after the initiation of transfer 
presented as a function of the initial charge Sq . 
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Fig. 11 — The fractional diminution of a one as a function of the amplitude of 
a "fat" zero in an alternating series of one's and zero's. This example also refers 
to a 25-fim 2-phase p-channel CCD. 
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the device, the left end follows, but the profile which is established 
relatively early in the device cycle remains at least superficially un- 
changed through the course of the transfer. 

If the potential of the source well is being changed at the same time 
the transfer is taking place, it is valuable to know the time develop- 
ment of the highest charge density in the source well, because it is 
necessary to prevent this point from ever reaching an injection con- 
dition. In this approximation injection occurs when Sq > y{V A — 
V F b) , and the well is no longer capable of holding the charge ; the well 
and the substrate then form a forward biased p-n junction. For the 
example discussed before, the 25-/im p-channel CCD, the highest volt- 
age in the well has been plotted as a function of time, assuming initial 
charges of 10 volts, 2 volts, and 1 volt; Fig. 13 shows the real time 
for the specified device and the general, nonnalized time. This figure 
shows that the potential on the plate must not drop to zero in a time 
less than 100 ns, if the injection of stored CCD charge into the bulk 
of the semiconductor is to be avoided. 

It is unlikely that the initial charge configuration in a CCD will 
be perfectly square as assumed. To test the effect of this assumption 
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Fig. 12— The charge density as a function of position under the CCD plate, 
and the parameter is the time after initiation of transfer. The percentages refer 
to the total fraction of charge remaining under the plate. 
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Fig. 13 — The amplitude of the largest surface potential difference 8q(0) as a 
function of time. Two time-axes have been provided, the generalized t, and 
/ for a 25-/j.m p-channcl device. 



on the results presented here, two initial conditions were used; one 
had charge distributed uniformly to an Sq value of 8 volts, and the 
other had charge distributed nonuniformly but having the same average 
charge density, that is Sq = 8 volts. Figure 14 shows the difference 
in the charge remaining in the potential well as a function of time for 
these two distributions, and it is easily seen that beyond the time of 
t = 0.25 the two charge quantities are substantially identical. 

The effect of finite carrier velocity at the transfer point was also 
tested. If infinite velocity cannot be achieved, the charge density at 
the point L cannot be 0, but it must be that finite value which will 
allow current to flow out of the potential well at the saturation 
velocity. Using this fact, it is possible to write the charge density 
at L as a function of the time rate of change of the total charge in 
the well dQ/dt. 



U (L) = Sq(L) = -- = --■ 



(21) 
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Here v 8 is a saturation velocity and, assuming its value to be 10 7 centi- 
meters per volt-second, the solution of equation (10) was reprogramed 
to reflect a value u(L) in accordance with equation (21). The dif- 
ference caused in the result was negligible, but because of the L 2 ///, 
time scaling, it was more significant for very small plates than for 
the larger plates. The time error is shown in Fig. 15 as a function of 
charge remaining beneath the plate. 

v. DISCUSSION 

We have developed an equation governing charge transfer in flat- 
plate CCD's and shown how it can be solved numerically by an un- 
conditionally stable finite difference routine. These solutions have 
been applied to the calculation of the properties of CCD's subject to 
certain constraints ; the most serious of these limitations is the assump- 
tion of the instantaneous creation of the transferring condition. This 
assumption transforms into physical reality as the use of square- 
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Fig. 14— Two initial conditions are compared to evaluate the impact of non- 
uniform distribution of a given total charge on the transfer predictions. The 
ordinate is the fractional difference in the total remaining charge and the 
abscissa is normalized time. 
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Fig. 15 — An illustration of the effect of limiting the carrier velocity at x = L 
to 10 7 cm/s, expressed as an error between the time predicted in the idealized 
case and the time in the velocity-limited case. 



wave driving. Under these conditions, we have shown that with 
empty zero's losses of the order of 1 percent can be achieved in the 
vicinity of 1 MHz in p-channel CCD's, and the same order of loss at 
5 MHz in n-channel CCD's. Further it has been shown that the use 
of "fat" zero's can effectively reduce this loss by as much as an order 
of magnitude. The problem which has not been solved here is that 
which arises when the accepting well is not created instantaneously but 
moves down with a finite time-rate of potential change. In this situa- 
tion, a one will have a longer time in which to transfer than a zero 
will. This may tend to reduce the amount of charge left behind by the 
one as contrasted with that left behind by a zero. This would diminish 
the loss. At this stage, it is evident that losses of realizable CCD's, i.e., 
those with plates in the range of 15 to 25 /xin, will be modest in the 
5 to 10 MHz range, particularly if they are made using an n-channel 
technology. Using properly loaded one's and zero's, short shift registers 
should deliver useful signals, even clocked at 50 MHz. 

This analysis has also shown that the existence of field-free regions 
in CCD's represents a considerable limitation on their performance. 
It will be possible to further lower losses by arranging substantial 
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penetration of transverse electric fields under the plates so that the 
charge transfer is limited by drift rather than diffusion. In this case, 
the scattering limited velocity assumes an important role. 

CCD's on 10 O-cm substrates appear, on the basis of preliminary 
measurements, 3 - 6 to behave according to the predictions of the space 
charge assisted diffusion theory. 
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{Added in prooj.) A related analysis with more limited scope than that above 
has been published by W. E. Engeler, J. J. Tiemann, and R. D. Baertsch in 
Applied Physics Letters, 17, 1970, p. 469. 



